An elastic, plastic, viscous model for slow shear of a liquid foam 
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We suggest a scalar model for deformation and flow of an amorphous material such as a foam 
or an emulsion. To describe elastic, plastic and viscous behaviours, we use three scalar variables: 
elastic deformation, plastic deformation rate and total deformation rate; and three material specific 
parameters: shear modulus, yield deformation and viscosity. We obtain equations valid for different 
types of deformations and flows slower than the relaxation rate towards mechanical equilibrium. In 
particular, they are valid both in transient or steady flow regimes, even at large elastic deformation. 
We discuss why viscosity can be relevant even in this slow shear (often called "quasi-static") limit. 
Predictions of the storage and loss moduli agree with the experimental literature, and explain with 
simple arguments the non-linear large amplitude trends. 



I. INTRODUCTION 

Elastic materials deform reversibly [l|; plastic materi- 
als can be sculpted, that is, they can be deformed into a 
new shape and keep it [|J ; and viscous materials flow Q . 
A wide variety of materials display a combination of these 
properties, such as elasto-plastic metals and rocks, visco- 
clastic polymer solutions or visco-plastic mineral suspen- 
sions 0,0113. 

Liquid foams, that is gas bubbles separated by liquid 
walls, are visco-elasto-plastic 0, [H, they are elastic 
at low strain, plastic at high strain and flow under high 
shear rate. This is also the case for other concentrated 
suspensions of deformablc objects in a liquid pi Hoi. fll| . 
such as droplets in emulsions, vesicles suspensions, or red 
blood cells in blood. 

Despite a large literature on experiments and simula- 
tions (see @ for a review) , we lack an unified theoretical 
description of foams. There is no consensus yet on a cen- 
tral question: what arc the physically relevant variables? 
A series of statistical models focus on fluctuations and 



Conversely, recent 
23l focus on average 



their correlations [121 1131 . 114. ,11 
contributions fl7l Iia Il9l |21|. |2 

macroscopic quantities to obtain a more classical contin- 
uous description. 

Here we choose to group three macroscopic quantities 
which are measurable as averages on microscopical de- 
tails [l?]]: (i) Elastic deformation is a state variable 24 1 
reversibly stored by the foam's microstructure, that is, 
the shape of bubbles [2^, [26J; it determines the elas- 
tic contribution to the stress, (ii) Plastic deformation 
results in energy dissipation analogous to solid friction, 
(iii) Large scale velocity gradients are associated with a 
viscous friction. Each of the three mechanical behaviors 
is associated with a material specific parameter: elastic 
modulus, yield deformation and viscosity. 

For simplicity, we assume here that these parameters 
are constant and the equations are linear. We consider 
here homogeneous deformation of a material, not depend- 
ing on space coordinates. We consider only the magni- 
tude of deformation, but not spatial orientation: the ma- 
terial state variables are all scalars. This represents an 



incompressible liquid foam, where the deformation is a 
pure shear. We assume that this shear is slow enough so 
that the foam is always close to mechanical equilibrium, 
but quick enough to neglect coarsening such as due to gas 
diffusion between bubbles, or bubble coalescence due to 
soap film breakage. Although this model is minimal, it is 
written with enough generality to enable for extensions 
to higher dimensions using tensors (the correspondancc 
with tensors introduces a factor 1/2, see section llV C[) . 
to higher shear rates, and to other ingredients such as 
external forces (to be published). 

This article is organised as follows. Section |TT] intro- 
duces a visco-elasto-plastic model (eqs. I3I7[) based on two 
scalar variables: the elastic deformation and the (slow) 
shear rate (Fig. [3]). The rate of plastic deformation is 
determined by both the applied shear rate, and the cur- 
rent state of the elastic deformation (or equivalently the 
elastic part of the stress) rather than by the total stress 
[27l . l2l| . Section IIIII presents scalar predictions of creep 
and oscillatory responses. The storage and loss mod- 
uli predicted as a function of the strain amplitude agree 
with experimental data without any adjustable parame- 
ters, using only the three model-independent parameters 
determined by experiments (yield point, shear modulus, 
viscosity). The agreement becomes very good if we de- 
scribe the plastic yielding as a gradual transition spread- 
ing between an onset value of deformation and a satura- 
tion value (eq. [5]). Section [IVI summarises and discusses 
our model, and opens some perspectives. 



II. MODEL 
A. Kinetics 

1. Elastic and plastic strain 

The elastic deformation U is a state variable, that is 
an intrinsic property of the foam's current deformation 
state. We note its time derivative dU/dt. Conversely, 
we use a dot for the total strain rate e and the plastic 
strain rate ip, emphasising that they are not the time 
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derivative of a state variable. For instance, the time in- 
tegral e — J e dt of the velocity gradient is the gradient 
of displacement (more generally, for large deformations, 
it is a function of the displacement): it is extrinsic and 
explicitly depends on the sample's past history. 

The total applied deformation rate is shared between 
clastic deformation U and the plastic deformation rate: 



dU 
dt 



E P . 



(1) 



In the particular case of an elastic regime, ip = 0, 
the elastic deformation U is equal to the total applied 
deformation on the material e. Thus, in an elastic regime, 
no intrinsic definition of U is necessary. 

However, as soon as ip ^ 0, the situation changes. 
U and e become independent variables, and e = J idt 
does not define the elastic deformation. In the extreme 
example of a steady flow, dll/dt = 0, then e = ip: U 
and i are no longer correlated. 

These variables are macroscopic: U is related to the 
clastic contribution to macroscopic stress and ip to the 
irreversibility of the stress versus total strain curve. In 
the specific case of foams, they can be traced back to 
detailed properties of the bubbles pattern: independent, 
intrinsic definition [24[ based on geometry (shape of bub- 
bles [2^ 1) for U- and topological rearrangements called 
"Tl processes" [lj, [29|, |3(| (using their rate and orienta- 
tion [17]) for ip. 



2. Sharing the total strain 

The problem now is to express how, in eq. ([Tj), i is 
shared between dU/dt and ip. We must write a closure 
relation between these variables, for instance by express- 
ing how ip depends on the current state of elastic de- 
formation and on the applied deformation rate: ip(U,i). 
We use the following three hypotheses leading to eq. O 

First we describe an abrupt transition from elastic to 
plastic regime, as could be the case for an ordered foam 
[311 ]. To indicate that Tls appear when the absolute value 
of deformation | U | exceeds the yield deformation Uy , we 
introduce the discontinuous Hcavisidc function TL (which 
is zero for negative numbers, and 1 for numbers greater 
than or equal to zero) . This hypothesis can be relaxed in 
the section III A3! introducing a more progressive transi- 
tion. 

Secondly, we account for the hysteresis. Plastic rear- 
rangements occur when the deformation rate e and the 
current deformation U have the same sign, and again we 
express it using TL. Else, the deformation rate results in 
clastic unloading, and the deformation gets smaller than 
the yield deformation. 

Thirdly, we use the fact that, in a slowly sheared mo- 
tion, the only relevant time scale to fix the rate of plastic 
rearrangements is i. 

Eventually, the plasticity equation writes: 



Eq. ([2]) can be used to close the system of equations. 
Injecting it in eq. flj yields an evolution equation of U 
as a function of the applied shear rate i: 



^-=i [l-TL(\U\-U Y )TL(Ui) 
dt 



(3) 



In eq. © Uy appears as the stable value for U, that 
is, a fixed point, at least if i > 0; else, the stable fixed 
point is —Uy- 



i P =TL{\U\-Uy)TL{Ui) i. 



(2) 



U 
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FIG. 1: Analog scalar system: an elastic brush whose flexion 
is U stick/slipping on wall. We represent several states, for 
an imposed oscillatory "painting-like" motion of the handle 
e, from rest position 0: (a) onset of sliding to the right, (b) 
far-right position, (c) far-left position. 

To visualise the direction and the amplitude of the de- 
formation U , we suggest an analogy with the motion of a 
brush on a wall (Fig. [TJ. The handle of the brush moves 
with an oscillatory position e parallel to the wall (ana- 
log of the imposed scalar deformation of the material), 
while the displacement of the handle with respect to the 
brush tip is U (the analog of the internal elasticity of 
the material). The sliding velocity of the contact point is 
therefore ip according to equation JTJ and is the analog 
of plasticity in a material. 



3. Gradual transition to plasticity 

In a disordered foam, for instance with a wide distri- 
bution of bubble sizes, topological rearrangements do not 
necessarily occur for the same value of deformation. 

We therefore distinguish two different yield deforma- 
tions. First, a plasticity yield U y , where deformation 
ceases to be reversible, as defined in material sciences. 
It is the highest deformation for which there is no Tl. It 
is characteristic of the microstructure, and can even be 
close to zero for a very disordered foam. 

Second, a saturation yield Uy, the saturation value of 
elastic deformation at which the material can flow with 
arbitrary large total deformations (for instance in Bing- 
ham fluids). It is the lowest deformation for which the 
Tls convert the whole total strain into plastic strain. 
That is, Uy is the collapse limit at which a material 
structure cannot sustain stress. 
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We interpolate between U y and Uy using a function 
h(U) which we call a yield function. It should be a 
growing (or at least non-decreasing) function of U for 
U > 0, and h(-U) = h{U). Moreover, h(0) = 0, so 
that h(U) > for all U. Beside that, there is no spe- 
cial requirement on h, which even does not need to be 
continuous. Now, U y is defined as the largest value of U 
for which h(U) = 0, and Uy as the smallest value of U 
for which h(U) = 1. They do not necessarily correspond 
to any singularity in h. We show in Appendix that the 
precise shape of h is unimportant: only U y and Uy de- 
termine material's behaviour. However it is useful for 
theory to derive analytical equations. 

The yield function h depends on the material under 
consideration, and can in principle be measured experi- 
mentally. By definition, h = corresponds to a purely 
clastic state where the elastic deformation follows the 
applied deformation. Conversely, for h = 1 the plastic- 
ity rate is equal to the deformation rate. Such a smooth 
transition from elasticity to plasticity generalises the pos- 
tulate ([2]) as: 



e P = h(U)H(Ui)s. 



(4) 



Note that we could in principle smoothen out the remain- 
ing Heaviside function too: depending on microscopical 
details, it could be conceivable that some Tls appear 
during the unloading. We do not explore this possibil- 
ity here, because we seldom observe this effect and it 
does not seem to improve significantly the predictions 
presented below. 

Injecting eq. (|4]) into cq. ([!]) we obtain: 



[l-h(U)H(Ui)\. 



(5) 



Again, the fixed points are U = ±Uy according to the 
sign of i. 

B. Dynamics 

1. Slow shear: foam close to equilibrium 

In a foam, bubbles can swap neighbours giving rise 
to Tl topological rearrangements. A Tl is an infinitely 
short event during which the energy is continuous. Thus 
it does not dissipate energy by itself, but it brings the 
foam in an out-of-equilibrium state. It is thus followed 
by a dissipation of energy during the relaxation towards 
another equilibrium configuration, with a microscopical 
relaxation time r ro i ax . 

The average life time of a contact between two bubbles 
is / , where / is the average frequency of Tls per bub- 
ble contact. If /Tjoiax *C 1, the foam evolves (it is not 
static) but spends most of the time at or very close to me- 
chanical equilibrium states. Thus Plateau rules of local 
mechanical equilibrium Q are (almost) always satisfied, 
up to corrections of order /r ro i ax . 



The frequency / can be determined by various causes 
of perturbations, for instance coarsening In rheology, 
it is determined by the plastic deformation rate ip pjj. 
For dimensional reasons, / is proportional to ip. Since 
the plasticity amplitude ep is always smaller than the 
deformation rate i (see eqs [5] and SJ) , the regime close to 
equilibrium is obtained in the slow shear limit: 



^Trelax < 1- 



2. Contributions to total stress 



(6) 



We now include an additional viscous dissipation from 
the global deformation of the network of bubbles. This 
contribution is not linked to the relaxation of rearrange- 
ments, and does not modify the slow evolution of defor- 
mation. 

We consider two separate contributions to stress, un- 
der the following hypotheses. According to experimen- 
tal tests [26|, we consider an elastic contribution to the 
stress <7 cl = [ill proportional to the elastic deformation 
U, where fi is the shear elastic modulus. It describes 
a classical elastic behaviour, with a reversible restoring 
force. 

According to the model proposed by Kraynik and co- 
workers @, [33 ] , we consider a viscous contribution to the 
stress due to large scale velocity gradients: ct v1s = r]i, 
where 77 is a macroscopic viscosity. It describes a classical 
fluid behaviour: the corresponding dissipated power is 
quadratic, proportional to e 2 . 
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FIG. 2: A linear elasto-visco-plastic Theological model. 

In the spirit of a polymeric model Q , we assume that 
the stresses add up (Fig. [2|): 



a = a + a 



= [iU + rji. 



(7) 



The material is characterised by the coefficients 77, (i 
and Uy (and U y in the case of gradual plasticity). Mea- 
suring experimentally, and understanding theoretically 
the physical origin of these coefficients, requires specific 
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FIG. 3: Scalar phase diagram for a slowly sheared foam or an 
emulsion. Axes are experimentally measurable [Th] local vari- 
ables: shear rate e and elastic deformation U. The crossover 
from elastic to plastic is defined as the onset of the first iso- 
lated topological rearrangements; it occurs around Uy, with 
possible precursors around U y . The yield deformation Uy cor- 
responds to a macroscopic rate of topological rearrangements. 
The crossover from solid to fluid is defined by the equality of 
viscous and elastic stresses. The slowly sheared regime pre- 
sented here ceases to be valid when e becomes comparable to 
r relax: inverse of the microscopical relaxation time. 



studies for each material considered: this is beyond the 
scope of the present paper. In principle, they can be 
rank- four tensors (anisotropic material). They can even 
vary with the material's state (non-linear material), for 
instance in a shear-thinning case. 

As opposed to the cross-over from elastic to plastic 
regimes, which is topological and is visible on images, 
here the crossover from the elastic to the fluid regime can 
be detected only by measuring forces. It occurs when the 
viscous contribution to the stress becomes larger than the 
clastic one. Fig. ([3]) thus plots the line corresponding to 
the crossover: [iU = r]i. 

Defining the macroscopic local Weissenberg number as: 



look the same; their difference is not kinematic but dy- 
namic: stresses arc independent on and proportional to 
i, respectively. 

3. Dissipation 

The close to equilibrium criterion (eq. [H]) regards time 
scales, and is not a criterion based on the absence or pres- 
ence of dissipation. Viscous dissipative effects can indeed 
be observed when considering measurements of the loss 
modulus at very low amplitude oscillations, and hence 
at very slow shear rate (as presented below in Figs. [9]- 
ITTT) . In fact, dissipation is absolutely necessary to relax 
towards equilibrium: it damps oscillations and decreases 
the energy. Note that a "quasi-static" regime, that is 
a succession of equilibrium states, necessarily obeys the 
equilibrium criterion; but it is not sure that the reverse is 
true. In fact, Ref. 33| claims that in the slowly sheared 
Couette flow by [34 1 the velocity profile is determined by 
the ratio of velocity-dependent forces (internal viscosity 
and external friction on the plates of glass): static simu- 
lations are inappropriate. 

As already mentioned, a Tl by itself, that is a side 
swapping, is an infinitesimally short topological event, 
during which the foam energy is continuous: there is 
no instantaneous dissipation. However, the Tl puts the 
foam in an out-of-equilibrium state. During a time r rc iax 
the foam relaxes to a local energy minimum by dissi- 
pating an energy 5E. A smaller microscopic dissipation 
yields a shorter relaxation time r re i ax , and a larger in- 
stantaneous dissipated power, of order 'Pdiss = 5E /t^i^. 
But the amount of energy dissipated, SE, is independent 
on the dissipation. Thus the dissipated power averaged 
over a long time (longer than e _1 ) is of order: 

(Pdiss) = fSE ~ ipSE ~ iSE. (9) 

This dissipated power is proportional to i, and not 
quadratic as in viscous flows, although the microscopi- 
cal origin is a local viscous dissipation [35| . 



Wi M = ^ 



(8) 



the cross-over between elastic and fluid regime occurs at: 
Wi M = CZ- 
Upon increasing e from the plastic regime where U is 
close to Uy , we predict a cross-over from a stress bounded 
by a constant value fJ-Uy (with a dissipated power linear 
in e, see next section) to a stress linear in r\e (with a 
dissipated power quadratic in e), characteristic of a vis- 
cous friction. The cross-over from a plastic regime to a 
fluid regime occurs when rye is equal to fiUy, i.e. when 
WiM = Uy. 

Since the plastic deformation is not bound, in the plas- 
tic regime the foam can flow indefinitely. As in hydro- 
dynamics, the displacement field itself is no longer rel- 
evant. The plastic flow 0, Q and the viscous flow Q 
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FIG. 4: Schematic impact of individual microscopic rear- 
rangements on the stored elastic deformation U , for a con- 
stant loading rate e. Rearrangements relax exponentially the 
deformation over a time r ro i ax , with here £T ro i ax = 0.02 <C 1. 
In the present macroscopic model, rearrangements are coarse- 
grained. 

The clastic deformation is almost independent on the 
shear rate e (see figure 2]). To obtain a steady shear in 
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a solid regime, when U saturates at the value Uy, an 
experimentalist has to apply an constant external force 
which balances the average elastic stress, and does not 
depend on e. 

Such a dissipated power linear in e, and a steady force 
which does not depend on e, are characteristic of a solid 
friction [MUzI- 



a = ay f(e/ey): the function / corresponds to the 
present elastic deformation U, while ey is the yield de- 
formation they chose equal to 1 for simplification. 



B. Creep under constant applied stress 



III. PREDICTION AND TESTS 

We model the foam response in one type of mechani- 
cal experiment, imposed deformation, and in two types 
of rheometrical experiments, creep flow and oscillating 
shear. 



A. Imposed shear 

Here, we calculate the transient response during a 
shearing experiment, that is, the relation U(e) between 
applied strain e and elastic deformation U. For simplicity 
we take here e = U = at the start of the experiment, 
but that assumption is easy to relax. 

By direct integration, see appendix, we show the mate- 
rial's response: the elastic deformation U is close to the 
imposed strain e at low applied strain, and tends to a 
saturation value at large applied strain. This robust be- 
haviour does not depend much on the chosen yield func- 
tion (see fig. [5}. 




e/U v 



FIG. 5: Response to imposed shear for two examples of yield 
functions. Dashed solid line: abrupt transition, h(U) = 



H(U - U Y ) (eq. 
terpolation h(U) 



0. Thick line: 
= (U-U V )/(U Y 



finite U y , 

~ Uy)H(U 



and linear in- 
-U y ) (eq.EJ. 



Here U y = 0.75 Uy- See fig. [12] for more examples. 

Thus the distribution of bubble sizes does not affect 
much the foam's transient response (as opposed to the 
liquid fraction, which drastically affects Uy [Uj]). This 
explains why in the literature the function U (e) is some- 
times taken for simplicity as a piece-wise linear function 
or as a hyperbolic tangent fill |. 

This provides both the physical origin for the func- 
tion (7(e) of the model by Janiaud et al. [HI], and a 
justification for their (up to now arbitrary) expression 




FIG. 6: Response to constant applied stress: (a) elastic defor- 
mation versus time and (b) local strain rate versus time, with 
a sharp transition to plasticity at U = Uy- In both figures 
the stress grows from bottom to top: below yield values, <7 a p P 
equals to 0.5, 0.7, 0.9 times ay (solid lines) and above yield 
values, cr a pp equals to 1.1, 1.3, 1.5 times \jUy (dashed lines). 
The time is adimensioned by A = rj/fj,. 



A creep experiment in a rheomcter applies a constant 
stress cr ap p . It determines the effective viscosity rj e f{ from 



the steady shear rate: 



VcS 



lim ^f. 

t^oo i(t) 



(10) 



The rhcological response is found from cq. ([7]) with a = 
cr ap p, and from eq. ([3]). The elastic loading and the strain 
rate are plotted on Fig. ([6]). The elastic deformation 
saturates to cr app //x when it is below the threshold Uy, 
and that it saturates to Uy when above the threshold, 
over a characteristic time A = rj/ '/j. 

At long times, the strain rate tends towards vanishing 
values below yield stress (the flow stops), and tends to fi- 
nite values above the yield: e(t — > oo) ~ (er a pp — fJ-Uy)/ rj. 




FIG. 7: Response to constant applied stress: (a) elastic de- 
formation versus time and (b) strain rate versus time, with 
a smooth appearance of plasticity in between a plasticity 
onset threshold U y and saturation Uy- Plasticity starts at 
U y = 0.75 Uy, with the linear interpolation of eq. 1191 Same 
legend as previous figure. 



We thus deduce the effective viscosity: 

r] cS = oo when cr app < [iU Y , (11) 

VeS = ~ 77/777 wncn O'app > IJ>Uy- (12) 

^app 

Taking a smooth plastic transition (eq. [5]) does not 
changes the overall features, except that the deceleration 
times below yield are no longer superimposed, see fig. 



C. Oscillating shear 

A Couette apparatus is another typical rheometry ex- 
periment 0, [ill, [2^, [3^] . It measures the stress a(t) on the 
walls while imposing an applied sinusoidal shear strain of 
pulsation cu = 2t:/T: 

e = j sin(u;t) . (13) 



-1 




FIG. 8: Long time periodic response to an oscillatory shear, 
for three amplitudes: e/Uy = 0.75, 1.1 and 1.8. (a) in the 
case where the plasticity appears abruptly at Uy , (b) gradual 
transition in between U y and Uy, with U v /Uy = 0.5 (eq. I19|l . 



1. Hysteresis cycle and non-linear response 

To test the effect of hysteresis of the model we calculate 
the response to such a periodic oscillatory shear bounded 
by amplitudes 7 and —7. In the slow shear limit, the 
frequency does not play any explicit role. We can thus 
keep it fixed without loss of generality. 

We integrate eq. ([3]). The periodic elastic deformation 
vs strain curve is plotted on Fig. ([H]) top. The stress 
response is linear in strain below the threshold, and sat- 
urates above in plastic regime, exhibiting a strong hys- 
teresis. 

Reversing the sign of the loading instantly stops any 
plasticity and the reponse becomes purely elastic. Multi- 
ple loading does not increase the slope of the loading part, 
nor the value of saturation yield; the foam is described 
as perfectly plastic. Such features are observed in exper- 
iments on other amorphous solids [1^, H(| (as opposed to 
strain- hardening features of crystalline metals [2|). 

Integration of eq. ^ describes a smooth variation of 
deformation, see Fig. © bottom. 



2. Storage and loss moduli: predictions 

In complex notation the stress response a* is linked to 
the strain e* by a* = (G' +iG")e* . Here G' is the storage 
modulus and G" the loss modulus of the material, defined 



7 



as the in-phase and out phase part of the response (first 
term in a Fourier series, see non linear models [l9 . l4ll . 

IH). 

When increasing the amplitude 7 of the imposed sinu- 
soidal shear strain, the response is first linear until the 
amplitude at which G' and G" start to vary. In both 



the linear and non- linear regimes, 
moduli arc calculated as: 



the storage and loss 



1 1 f 

G = ? — / a{t)di, 

7 ^ Jo 
1 1 f T 
G = / a(t)de, 

7 * Jo 



(14) 



G' is proportional to the area enclosed by the (cr(i), s(i)) 
curve, while G" is proportional to the area enclosed by 
the (a(t),e(t)) curve. When plasticity occurs, the cy- 
cle has a non-vanishing area in the (cr(t),e(t)) diagram, 
meaning a non- vanishing loss modulus G" . 

In the present model a(t) depends on the current elas- 
tic deformation U(t) and shear rate e(t) (eq. \7§. For 
the case of the abrupt elastic/plastic transition, the an- 
alytical integration of areas is simple and provides the 
following solutions for the moduli. Using cqs. (| 14[) we 
obtain, when 7 <C Uy- 



G' ~ n 
G" = 77, 



(15) 



which is the usual linear visco-elastic regime. Note that 
our model predicts frequency independant moduli, for a 
fixed small amplitude 7. At large amplitudes, when 

7 > U Y : 



3/2 



G' ~ n- 

77 V 7 J 

AU Y 

n — +v- 

7T7 



G" 



(16) 



These asymptotic dependencies in 7~ 3 / 2 and 7 _1 , re- 
spectively, are obtained analytically and are robust with 
respect to the model. The analytical expression of G' and 
G" over the whole range of 7, but with 77 = is explicitly 
presented in refs. [2ll. |22|. 

For a smooth yield function h, predictions are obtained 
numerically and plotted on Figs. ([9l fTT|) . 



3. Comparison with experiments on emulsions and foams 

Rheometry measurements of monodisperse emulsions 
[H| (Fig. El) and polydisperse foams [33, [H| (Figs. \M 
and [TTj) directly yield, without hypotheses, the values 
of the material parameters required by the model. The 
shear modulus /1 is read from the value of G' at low am- 
plitude. The viscosity 77 is read from the value of G" at 
low amplitude (or the value of the minimum, in Fig. [9l 
where the two data points at 7 < 10~ 3 have too large 
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FIG. 9: Storage and loss moduli versus strain amplitude for 
a monodisperse emulsion. Symbols: experimental G' (circles) 
and G" (triangles) in a close-packed emulsion (Fig. 1 of ref. 
[43l ]. fraction of the continuous phase 20%, droplet size 0.53 
/im, oscillation pulsation w = 1 rad s" 1 ). Lines: models for 
G' (solid line), and for G" with an abrupt transition (dashed 
line), with viscosity (dash-dotted line), with viscosity and a 
smooth yield function h — (U/Uy) 2 , U y — (solid line). 
Model parameters: shear modulus fx = 1.7 10 3 Pa, yield de- 
formation Uy — 0.045, viscosity 77 = 30 Pa.s. 




FIG. 10: Same as Fig. Q for a polydisperse foam [4J]. Liq- 
uid fraction 5%, bubble size 40 to 70 /im, u — 1 rad s _1 . 
Lines: models for G' (solid line), and for G" with an abrupt 
transition (dashed line), with viscosity and a smooth yield 
function h = {U/Uy) 2 , U y = (solid line). Model parame- 
ters: fj, = 100 Pa, U Y = 0.2, 77 = 15 Pa.s. 
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FIG. 11: Same as Fig. for a monodisperse foam [3a |. 
Liquid fraction 8%, bubble size 21 fim, u — 1 rad s~ . Data 
of G' and G" are normalised by /i and 7 by while rjui/^i = 
0.08. Same legend as figure [lOl Inset: Strain-rate frequency 
superposition: same data plotted as a function of pulsation 
u>, for a given maximum strain rate 70 = to-y [20fl . 



error bars to be taken into account, according to T. Ma- 
son, private communication). The yield deformation Uy 
is read from the intersection of low amplitude plateau of 
G' and its large amplitude —3/2 exponent power-law. 

A purely elasto-plastic model is enough to predict G' 
correctly, over the whole range of amplitude, including 
the —3/2 exponent power-law. This simplest model also 
describes correctly the large amplitude trend for G" . 

The low amplitude value of G" can be modelled by in- 
cluding a viscosity (dash-dotted line of figure [5]), which 
confirms that viscosity is relevant even in such slowly 
sheared models. This procedure is at the expense of a 
slight over-prediction at large amplitudes. This latter 
aspect suggests a possible shear-thinning, that is a de- 
crease of the viscosity 77 with the shear rate, similarly to 
the observed reduction of the drag of foams in motion in 
channels [3l| . 

The agreement between the data and the model, with- 
out adjustable parameter, is good. It is still improved, 
even for G" at intermediate amplitudes, if we account for 
the fact that the first Tls appear gradually at a value U y 
lower than Uy (solid lines on figures for G"). Even the 
value of U v itself is not very important, and in order to 
avoid introducing a free parameter we use here U y = 
and a smooth yield function. 

These result suggest that data can be rescaled with the 
yield deformation Uy, and it suggest a rescaling when 
plotting data as a function of frequency, followin g th e 
strain-rate frequency superposition method (SRFS [2CJ|). 
This method consider measurements with a fixed max- 



imum strain rate 70 = luj, it is therefore equivalent 
to vary frequency or oscillation amplitude. The nat- 
ural rescaling for pulsation that appears is tu/b^o) = 
w/(7o/LV) = Uy/j, see inset of Fig (fTTj) , and our model 
predicts the global shapes of the moduli curves as ob- 
served in (20| . The characteristic frequency 6(70) is here 
linear in 70, also in agreement with the trend observed 
in [20( for large enough strain rates. 



IV. DISCUSSION 
A. Discussion of the predictions 

Independently from us, Hohler et al. solve a purely 
elasto-plastic model [2l], [22| ■ Since they neglect the vis- 
cosity, they can eliminate U and replace it by cr//i. The 
agreement of their model with data of Fig. [Til as we U as 
other experimental data, is good for G' over the whole 
range of amplitude; and also for G" but only at large am- 
plitude, where the dissipation comes from the relaxation 
after Tls rather than from 77. 

Our models describes better: (i) G" at low amplitude 
using viscosity; and (ii) G" at intermediate amplitudes, 
if we account for the fact that the first Tls appear at a 



v 
con- 



value Uy lower than Uy. 

The predicted curves are robust with respect to U 
This implies that we do not need to fit it; but that 
versely, we are not yet able to deduce U y from G" data. 
If we had a direct experimental measurement of U y , we 
could inject it in the model to predict G" at intermediate 
amplitudes, near Uy, but the resulting predictions would 
be very similar to the present ones. 

While U y corresponds to the onset of isolated plastic 
events, at a deformation Uy the plastic events have a 
macroscopic effect: they catch the total strain; the foam 
flows without increasing its deformation any longer. 

The effective viscosity rj e f[ diverges to infinite values 
when (T ap p//i — > Uy . This means that the foam comes 
close to its yield deformation, in the fluid sense. This 
change in behaviour was shown by [4a |. and modelled 
by a granular model with a velocity dependant friction. 
Here, the dynamics is entirely driven by a constitutive 
equation with shear-rate independent parameters. It is 
the transient clastic loading that drives a transient flow, 
which stops when the stress is not strong enough. 

The trends of figure ([6]) agree with experimental data 
on various material reported in (45j , at least with granu- 
lar materials and emulsion. The present model does not 
predict the apparent shear-thinning behaviour observed 



45 



where an 
The low 



with their experimental data with foams 
increase of shear rate with time is found 
amplitude predictions for the visco-elastic regime (see Eq. 
I15[) are frequency independant moduli. The limitation of 
the present model is thus that it does not fully capture 
the slight increase of the loss moduli at low frequency, 
and the uj 1 / 2 trend at large frequencies (see for instance 
the case of foams 0). 
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B. Elastic, plastic, viscous model 



A complete model for an elastic, plastic, viscous foam 
requires to recognise the role of three physical variables 
U, i, ip. There is a relation between them (eq. H}. 
Unless specific approximations apply, a foam's represen- 
tative volume element (RVE) is characterised by two in- 
dependent variables: we suggest to select the local clastic 
deformation U, and the local shear rate e, which are in- 
tuitive and physically relevant. Both of course depends 
on the sample's past history, but this history plays no 
explicit role. Both are always defined, whether in elastic, 
plastic or viscous regime [25[ • Two recent works [2(], E3] 
find that G' and G" depend on the strain amplitude and 
on the strain rate (rather than on the frequency), in the 
same spirit as our phase diagram (fig. [3]). 

Each volume element can thus be plotted as a point in a 
phase diagram (Fig. [3]); that is, the (e, U) plane [la ]. In a 
heterogeneous flow, different volume elements of the same 
foam are plotted as different points. A volume element's 
evolution is a trajectory on this plane. Simple materials 
correspond to the axes of the plane: pure elastic and pure 
plastic regimes on the vertical axis, pure viscous regime 
(Navier-Stokes) on the horizontal axis. 

The importance of U is the most original feature of 
the present model: U cannot be entirely determined by 
i since the latter can change sign; U cannot be entirely 
determined by a if the viscous contribution is not negli- 
gible. 

The yield function describing the occurrence of plas- 
ticity can be linked to the traditional hardening modu- 
lus, used for the description of plastic materials [48[. It 
is defined as K = da/de p , while the elastic modulus is 
fx = da /U . In the present model, the hardening modulus 
is dependant on the elastic deformation: K = 1/ h(U) — 1. 
It therefore vanishes when U tends to its saturation value 
Uy- at this point the material does not harden any more. 

A deformation beyond Uy is not accessible when start- 
ing from rest (Fig. [3]). But the foam could initially be 
prepared (for instance artificially [49j ) in a configuration 
very far from equilibrium. Under a steady shear rate e, 
the deformation U always tends towards Uy(s), whether 
from below or from above. 

In a flowing foam, there is always a viscous dissipation. 
Its contribution becomes dominant in front of the plastic 
dissipation if i > fiU/rj. This is compatible with the slow 
shear criterion, i <C (ie on a second Weissenberg 

number that is here Wi TO = £T rc i ax <C 1), if there is a 
scale separation between the microscopic relaxation time 
TYciax towards local equilibrium and the large scale time 
i]/fi. The dimcnsionlcss ratio yUT re iax/?7 of microscopic 
to macroscopic times is analogous to the parameter / of 
granular materials [50j . 



C. Perspectives 

This model could in principle be generalised to higher 
velocity gradients [23[. This would require a high flow 
velocity varying over a small scale, and r ro i ax could play 
an explicit role. The deviation from equilibrium, of order 
ETreiax, would become significant: for instance, under a 
steady shear the limit value of U could become larger 
than Uy- 

A rheometer such as a Coucttc apparatus can mea- 
sure (Ti2 = a (tangential force per unit wall surface) and 
E\2 = e/2 (components of the symmetrized deformation 
gradient), in a coordinate system aligned with walls. For 
comparison with tensorial data it is especially important 
to bear in mind that there is a factor 1/2: the threshold 
Uy on oscillation amplitude e as measured by a Coucttc 
rheometer, corresponds to a threshold Uy/2 on the ten- 
sorial deformation £12. The present scalar approach can 
be generalised to take into account such an influence of 
the orientation of material deformation, as well as spatial 
variations [13, [EH . 

The present paper is a contribution to a lively debate. 
Can a foam be described as a continuous material? We 
tend here to answer "yes" , in the same spirit as many 
recent papers which describe or predict rheological prop- 
erties at large scale [H [U, HH, [H [H H, [E| . Statis- 
tical de scrip tions of fluctuations and their correlations 
[I2T [HI, Il5l . [iH [45| are then useful in describing even 
more complex rheological behaviour such as shear band- 
ing [3| or growing correlation lengths cale near the glass 
transition [531 ] - Interestingly, even in granular materials, 
where very large scale fluctuations are known to occur, a 
recent paper emphasises the dominant role of the "con- 
tinuous material" description based on averages [50l ] . 
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APPENDIX: ELASTIC-PLASTIC TRANSITION 

1. Transient response from rest 

We assume (in this section only) that the deformation 
rate e keeps a constant sign. Under this essential assump- 
tion, we can calculate analytically the transient response 
during a shearing experiment. That is, the relation U(e) 
between applied strain e = J e dt and clastic deformation 
U. 
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The yield function h is defined to interpolate between 



h(U) = for 
integration, eq. 



< U < U y , and h(Uy 



yields: 



dU 



1. By direct 



(17) 



l-h(U)' 

Here, without further loss of generality, we have also as- 
sumed (but it is easy to relax) that e = U = at the 
start of the experiment, and that e > 0, so that U > 
too. 

Eq. p7|) yields the function e(U), which can be in- 
verted to obtain U{e). These functions can be measured 
on experiments and compared with predictions derived 
from direct measurements of h(U). 

Whatever the function h(U), eq. (fT7|) implies that U w 
e as long as U < U y : applied and elastic deformation are 
equal in the elastic regime. At the onset of plasticity (or 
topological changes), U > U y , they differ. When U gets 
close to Uy the r.h.s. of eq. fTTj) diverges. Thus, when e 
increases arbitrarily, U tends asymptotically towards the 
saturation value Uy. 

2. Examples of yield functions h 



Yield function h Elastic response U(e)/Uy 



H(U - Uy) 

finite U y (eq. [T9j 

(U/U Y )° 

(U/Uy) 1 

(U/Uy) 2 

sin 2 {U/Uy) 



e - U{e/U Y 
Eq. J20J 


1 — exp(— e/Uy) 

tanh(e/t/y) 

arctan(e/C/y) 



TABLE I: Elastic deformation for different examples of yield 
function h, for a non-deformed initial condition U(0) = and 
with e of constant sign. 

Table [J proposes a few examples of yield functions h, 
and some are plotted on figure [T2"l 

Eq. ^ is only a particular case of the more general eq. 
([5]). with h being the discontinuous Hcaviside function: 



h(U)=H(U -Uy). 



(18) 



Eq. (fT7|) thus includes the case of the abrupt transition. 

An example of a yield function with finite U y is a piece- 
wise linear function: 



U<Uy 
U > Uy 



h(U) 



h(U) = 0, 
u-u v 



U Y -Uy 



(19) 



and eq. (jTTJ) yields directly: 



£ <Uy 

e>U v 



U(e) 
U(e) 



U Y -U y 
Uy 

UY-Uy 



( '-"v \ 

3 \U Y -Uy) 



(20) 




u/u, 



(a) 




e/U, 



(b) 



FIG. 12: Responses from rest for some examples of yield func- 
tions in Table I. (a) h(U); (b) U/Uy versus e/Uy- since they 
are very similar, for clarity only some of them are plotted. 
Thick solid line: abrupt transition, h(U) = Tt(U — Uy). Thin 
line: finite U y , here U y = 0.75 Uy, and linear interpolation 
h(U) = (U -U y )/(U Y -U y )H{U - U y ) (eq. fJU). Dashes: 
vanishing U y and quadratic interpolation, h = (U/Uy) 2 - 
Dots: vanishing U y and linear interpolation, h = (U/Uy)- 



We can interpolate between abrupt and smooth tran- 
sitions, using the family of model power-law yield func- 
tions: 



h(U) 



u 

~Uy 



(21) 



For instance, the quadratic expression h(U) = (U/Uy) 2 
yields U(e) = Uyt&nh(e/Uy). With these functions, 
plasticity appears more or less gradually, as soon as 



U > 0. That is, U v = 0. The limit 



oo is the 



Heaviside function (eq. [T8|) . 

More generally, the yield function can be thought as 
the convolution of the Hcaviside function TL and a distri- 
bution of yield values py : 



h(U) = / PY (Uy)H(\U\ - Uy)dUy 



(22) 



For instance, if the distribution of yield values p is a Dirac 
peak at Uy, it results in a Heaviside yield function h (eq. 
HE]). 



3. Robustness with respect to the choice of h 

Some functions U(e) from table U are plotted on figure 
112b . Strikingly, they do not depend much on the actual 
expression of h(U). In fact, only the expression of h near 
Uy matters; the relation between e and U is robust. The 
clastic deformation U is close to the imposed strain e at 
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low applied strain, and tends to a saturation value at 
large applied strain. 

The only important feature of h is its derivative h! just 
below the yield point: 

It determines how the fraction in the r.h.s. of eq. (JT7J) 



diverges. Thus U (e) is not the same if h! is zero or infi- 
nite, or even not defined as in eq. I|18p . If it is infinite, U 
reaches the saturation value at a finite value of applied 
deformation. 

Conversely, if h' is finite, as in most examples of Table 
I, the behaviour is universal. In eq. (fTT)) . the fraction 
diverges as (Uy ~ U)h' . Thus, whatever the value of h', 
e(U) diverges logarithmically and U(e) tends exponen- 
tially towards Uy. 
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